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ABSTRACT 

Semiconductor  alloys  like  Hgi_^CdyTe  and  Ali_^Ga^As,  where  there  is  a  close  lattice  constant 
match  between  the  constituents,  are  nearly  random.  However,  for  mesoscopic  size  scales  of  radius  r,  that 
are  large  compared  with  a  lattice  constant  (~25  A  <  r  <~  250  A),  hut  small  compared  with  typical  device 
dimensions,  the  number  of  substitutional  sites  is  small  enough  so  the  root-mean-square  concentration, 
(Ax)^  is  sufficiently  large  to  produce  random  arrays  of  mesoscopic  “quantum  dots”  that  can  adversely 
affect  many  device  properties.  The  influence  of  the  adverse  effects  differ  among  various  properties — for 
example,  electron  and  hole  mobilities,  lifetimes,  and  so  on — but  in  general  are  expected  to  be  worse  the 
smaller  the  fundamental  gap  becomes.  These  kinds  of  fluctuations  are  suppressed  in  lattice -mismatched 
alloys  like  //gi.x  ZnJ'e  and  Gai.x  In,As  because  there  is  a  long-range  strain  energy  penalty  associated  with 
them.  Thus,  lattice-mismatched  alloys  are  more  spatially  uniform  than  lattice -matched  alloys.  In  this 
paper,  a  calculation  of  the  impact  of  concentration  fluctuations,  among  regions  large  enough  so  it  is 
appropriate  to  assign  them  a  band  structure  (-390x780x780  A^),  on  tunneling  currents  is  presented.  An 
enhancement  of  -1.20  (-20%)  in  the  average  tunneling  current  density  for  LWIR  material  is  found,  and 
its  root  mean  square  deviation  is  -0.90.  For  25  pm  sized  pixels  this  translates  to  0.3%  variations  among 
pixels  due  to  this  effect.  This  estimate  does  not  include  the  impact  of  smaller  regions  that  must  be  treated 
quantum  mechanically.  As  pixel  sizes  decrease  the  impact  of  this  variation  will  become  more  significant. 
Surprisingly,  the  effect  of  fluctuations  on  the  leakage  current  characteristics  of  VLWIR  (20  pm  cut  off) 
pixels  is  also  small. 


Alloys  of  the  semimetal  HgTe  with  the  semiconductor  CdTe  and  with  other  II- VI  semiconductors 
can  have  their  band  gap  engineered  by  varying  the  proportion  of  the  constituent  compounds.  For  many 
applications,  such  as  the  use  of  such  alloys  as  detectors  in  infrared  focal  plane  arrays,  it  is  important  that 
the  devices  of  the  array  (pixels)  be  uniform,  and  that  the  material  within  each  pixel  be  homogeneous. 
Variations  in  alloy  composition  within  a  pixel,  and  among  pixels,  can  degrade  system  performance. 

The  lattice  parameters  of  ZnTe  and  CdTe  are  significantly  different,  unlike  HgTe  and  CdTe, 
which  are  almost  perfectly  lattice  matched.  The  energetic  consequences  of  this  difference  tend  to  reduce 
composition  fluctuations  in  the  Zn-bearing  alloy. 

The  lattice  parameters  of  most  ternary  III-V  or  II- VI  semiconductor  alloys  of  composition 
Ai_xSxG  follow  Vegard’s  law  quite  accurately,  with  little  bowing.'  When  the  compounds  AC  and  BC  are 
lattice  matched — for  example,  Al\_Aja,^As  or  Hgi_^CdxTe — the  cations  A  and  B  are  randomly  distributed 
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on  the  fee  cation  sublattice,  and  deviations  of  the  local  composition  from  the  average  x  are  determined 
purely  by  statistics.  When  there  is  a  lattice  mismatch,  any  region  of  the  alloy  whose  composition  deviates 
from  X  would  have  a  corresponding  unstressed  lattice  parameter,  and  therefore  must  be  either  compressed 
or  dilated  in  order  to  fit  into  the  space  available  to  it.  In  a  paper  which  linked  the  mixing  enthalpy 
parameter  of  such  alloys  with  their  bulk  elastic  properties, it  was  suggested  that  the  associated  elastic 
energy  might  reduce  composition  fluctuations,  and  the  suggestion  was  quantified  for  A/i_xGaxA5  and 
Ini^ypaJ^s.  In  a  later  paper,  the  magnitude  and  density  of  concentration  fluctuations  in  Hgi^J^dyTe  and 
were  presented."* 


In  the  absence  of  any  elastic  strain,  the  mean  square  of  a  deviation  Ax  from  the  average 
composition  x  of  a  region  containing  N  cations  ions  is 

=  . 

N 

The  elastic  energy  density  of  such  a  region  is^’^ 
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where  a  and  b  are  the  lattice  parameters,  and  B  is  the  bulk  modulus  of  the  alloy. 
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This  energy  and  the  mixing  entropy  of  the  alloy  were  used^  in  the  thermodynamic  formula  for  the 
macroscopic  composition  fluctuations  to  show  that  the  mean  square  fluctuations  in  a  lattice -mismatched 
alloy  are  reduced  below  those  in  a  lattice -matched  alloy  by  a  factor 
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where  V  is  the  molar  volume  of  the  alloy,  R  is  the  gas  constant,  and  T  is  the  temperature.  The  appropriate 
value  of  T  is  the  growth  or  annealing  temperature  at  which  the  alloy  reaches  equilibrium  before  cooling. 


In  the  application  we  are  considering  here,  the  primary  concern  is  the  occurrence  of  defects  which 
modify  the  performance  of  a  pixel.  The  defects  discussed  here  are  not  defects  in  the  usual  sense,  such  as 
impurities  or  vacancies,  but  regions  in  the  absorption  layer  or  junction  plane  whose  composition  (and 
therefore  whose  band  gap)  differs  from  the  average. 


In  a  random  alloy,  the  distribution  of  cation  populations,  P(x),  is  binomial;  the  binomial 
distribution  is  closely  approximated  by  a  Gaussian  distribution  for  as  few  ions  as  the  12  ions  comprising  a 
single  nearest  neighbor  cation  shell. 
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where  Ax  =  x  . 


Figure  1  shows,  as  deduced  in  Reference  4,  the  average  distance  between  defects  as  a  function  of 
AXniin,  the  departure  of  the  composition  from  the  average,  for  defects  containing  216,  1,730,  5,830, 

13,800,  and  27,000  cations,  in  Hgi_^CdjcTe  with  x  =  0.225  and  in  Hgi_^ZnJ"e  with  x  =  0.18.  For  regions  of 
Hga,22Cdo.nTe  with  a  size  (diameter)  of  ~  60  A  and  Ax^n  >  1.7%,  the  separation  between  the  mesoscopic 
regions  exceeds  their  size,  so  they  do  not  overlap  and  can  be  treated  as  isolated.  Smaller  regions  or  lower 
AXniin  overlap  and  are  meaningless.  While  the  fluctuations  of  regions  with  size  scales  >  60  A  are  real, 
these  regions  cannot  be  treated  as  though  they  have  a  distinct  band  gap  unless  their  size  exceeds  a 
de  Broglie  wavelength  )idB- 
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The  energy  we  have  used  to  specify  )idB  is  the  band  gap  Eg.  For  regions  smaller  than  PidB  a  full  quantum 
mechanical  treatment  is  required.  Let  us  concentrate  for  now  on  size  scales  larger  than  )idB- 


The  immediate  concern  is  on  tunneling  leakage  currents.  Other  properties  relevant  to  device 
performance  are  impacted  by  fluctuations  but  will  not  be  treated  here.  The  tunneling  rate,  T(x)  across  a 
reverse  biased  junction  is:^ 
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where  m*  is  the  electron  effective  mass,  e  is  the  electron  charge,  and  the  relevant  field  (E)  is  the  band  gap 
divided  by  the  depletion  layer  thickness  L, 

E  =  Eg/L  .  (7) 


If  the  junction  interface  is  divided  into  regions  that  are  one  a  de  Broglie  wavelength  thick  and 
(2/idB)^  in  cross  sectional  area,  then  the  average  tunneling  rate  across  these  regions  is. 


where,  for  lattice  constant  a. 
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Figure  1.  Concentration  Fluctuations,  Xc  -  10  |J,m. 


(a)  Hgi.^CdJ'e,  x  =  0.26.  The  average  separation  D[A]  between  spherical  regions  of  radius  r[A]  with 
magnitudes  greater  than  a  specified  minimum  x^in,  plotted  against  x^n- 

(b)  Hgy.yZnJ'e,  x  =  0.18,  equilibrated  at  T  =  220°  C.  The  average  separation  D[A]  between  spherical 
regions  of  radius  r[A]  with  magnitudes  greater  than  a  specified  minimum  x^n,  plotted  against  x^n. 


The  factor  of  1 6  arises  because  there  are  4  cations  per  unit  cell  and  the  volume  of  the  mesoscopic 
regions  being  used  is  4  The  gap  energy  as  a  function  of  x  can  be  approximated  by 
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Inserting  Eq.  (8)  into  Eq.  (11)  and  integrating  yields, 
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where  a  is  defined  as, 
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The  root  mean  square  tunneling  rate  can  also  he  found  in  closed  form,  and  is, 

T{x)  ^  ’ 


(14) 


Einally  the  current  I,  through  a  pixel  of  size  w,  is  the  sum  of  the  currents  through  the  array  of 
mesoscopic  regions  within  the  pixel,  so  we  have,  (I  /I(x)  =  (T  /T(x).  The  deviations  in  the  currents 
through  each  mesoscopic  region  add  in  quadrature  so  ((AI)^  =  (the  number  of  mesoscopic  regions  per 
pixel)  ’^^  ((AT)^  ,  which  is, 

{(A/)-)‘"  ,  ((aD-)- 

3 - — = - = - 3 - L —  (15) 

I{x)  {wUXJ  r(T) 

As  long  as  w  »  /idB  {-  ~390  A  for  10  |am  cutoff  material),  the  variation  in  the  leakage  currents 
from  pixel  to  pixel  will  he  modest. 


DISCUSSION 

The  results  for  various  quantities  are  collected  in  Table  1.  Aside  from  the  numbers  reported  in 
Table  1,  the  quoted  results  depend  on  the  following  parameters,  dEg/dx  =  1.7[eV],  L  =  l[(tm],  and 
a  =  6.5  A.®  If  the  carrier  concentration  or  the  operating  temperature  causes  L  to  change  then  the  results 
will  change  appreciably  since  a  is  proportional  to  l}  and  a  enters  in  an  exponential .  Eor  the  numbers 
used,  the  effect  of  concentration  fluctuations  on  the  average  leakage  current  for  -  10  pm  is  only  a 
20%  increase.  Surprisingly,  the  increase  in  the  leakage  current  for  Xco  =  20  pm  is  smaller,  7%.  This  occurs 
because  the  de  Broglie  wavelength  for  the  20  pm  cutoff  material  is  larger  causing  the  size  of  the 
mesoscopic  regions  considered  to  be  larger.  The  smaller  sized  regions  still  contribute  to  the  leakage 
current  even  though  this  analysis  cannot  properly  treat  their  influence.  To  estimate  their  effects  the  size  of 
the  regions  was  decreased  by  a  factor  p  in  the  analysis.  We  expect  that  this  over  estimates  the  impact  of 
these  smaller  regions,  since  quantum  effects  are  expected  to  smooth  the  band  structures  of  these  regions. 
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The  size  reduction  results  are  collected  in  Table  2.  Notice  that  as  r]  decreases,  ((AI)^  first  decreases, 
and  then  begins  to  increase.  The  initial  decrease  occurs  because  (w/2?Icib)  in  the  denominator  of  Eq.  (15) 
increases  faster  than  ((AT)  .  For  smaller  T],  this  situation  reverses,  and  ((AI)  increases,  although  it 

never  gets  very  big. 

Table  1.  Quantities  that  Enter  the  Fluctuation  Calculations  for  Two  Cutoff  Wavelengths 


Xco(lim) 

X 

Eg[mc  v] 

m* 

?ldB[A] 

N[10*’] 

E[V/cm] 

(r> 

T(x) 

/  o\l/2 

W) 

T(x) 

((A/r)‘" 

lix) 

10 

0.22 

0.12 

0.008 

390 

35 

1200 

1.20 

0.90 

0.0056 

20 

0.19 

0.62 

0.006 

630 

15 

620 

1.07 

0.41 

0.0032 

Table  2. 


Effect  on 


/  I{x)  of  Decreasing  the  Size  of  the  Mesoscopic  Regions  by  a  Factor  T] 


11 

Xco  =  10  p-m 

20  p-m 

1 

0.56 

0.32 

1/2 

0.45 

0.25 

1/5 

1.50 

0.41 

1/10 

12.00 

0.78 

1/20 

— 

2.80 

Fluctuations  also  have  consequences  on  the  bulk  properties  of  these  alloys.  For  example,  the 
r  =  3 1 .6  A  regions  with  Axmin  =  1 .7((Ax)^  constitute  a  random  array  of  “quantum  dots”  in  the 
conduction  band  of  HgCdTe  10  |am  cutoff  material  with  amplitudes  of  ±20  meV  and  average  separations 
of  ~70  A.  The  Axmin  =1.7((Ax)^  value  was  chosen  as  an  example  because  the  average  separation 
between  such  regions  is  larger  than  the  size  of  the  region,  so  there  is  no  issue  about  overlapping  regions. 
For  a  carrier  concentration  of  ~5  x  10^"^  cm“^  in  such  a  rough  potential  arrangement  at  0°  K  (assuming  the 
donor  states  resonate  in  the  conduction  band  so  there  is  no  carrier  freeze  out),  states  up  to  ~  30  meV  will 
be  filled.  This  estimate  assumes  that  half  the  volume  are  valleys,  the  electrons  fill  only  the  valley  states, 
and  takes  some  account  of  the  nonparabolic  nature  of  the  conduction  band  structure.  At  77°  K  the 
electrons  will  be  spread  over  an  additional  ~7  meV.  To  propagate  in  this  mountain  range,  electrons  must 
either  partially  tunnel  through  or  percolate  around  the  hills,  or  at  very  low  temperature  experience 
“Anderson  localization”  and  “Mott  variable  range  hopping.”’  In  any  case,  their  mobility  will  be  smaller 
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than  if  they  experience  normal  drift  transport.  This  behavior  in  low  n-type  HgCdTe  is  well  documented,* 
hut  has  never  been  ascribed  to  concentration  fluctuations.  From  Figure  1(b)  it  can  be  seen  that  the  average 
separation  between  fluctuations  in  HgZnTe  of  magnitude  1.7  %  is  ~1500  A  rather  than  the  ~70  A  for 
HgCdTe  with  the  same  band  gap.  If  the  HgZnTe  is  equilibrated  at  100°  C  rather  than  220°  C,  then  the 
separation  increases  still  further  to  ~1  pm. 

When  material  has  a  mobility  edge,  where  states  above  a  critical  energy  have  a  high  mobility  and 
those  below  this  energy  have  a  low  mobility,  it  leads  to  1/f  noise.®’^  Thus,  the  concentration  fluctuations 
are  also  a  likely  source  of  some  of  the  observed,  but  unexplained  1/f  noise  found  in  systems  made  from 
these  materials. 

The  fluctuations  will  impact  many  other  properties  important  to  device  applications.  The  list 
includes  hole  mobilities,  Shockley-Read-Hall  lifetimes  (through  capture  cross  sections  and  ladder  states 
in  the  gap),  and  Auger  lifetimes  (estimated  reductions  of  factors  of  3  at  higher  temperatures  to  10  at  lower 
temperatures  have  been  made).'^ 

While  HgZnTe  alloys  have  been  grown  by  liquid  phase  epitaxy'^  and  molecular  beam  epitaxy^* 
methods,  and  some  properties  characterized,  no  measurements  have  been  made  that  have  a  direct  bearing 
on  the  predictions  of  this  paper. 

Fluctuations  can  degrade  the  properties  of  lattice -matched  alloys  more  than  those  of  lattice- 
mismatched  alloys  where  they  are  suppressed.  The  magnitudes  of  the  effects  on  various  properties  must 
still  be  computed,  but  it  is  clear  from  the  qualitative  arguments  made  here  that  they  will  be  detrimental. 
Their  suppression  in  HgZnTe  should  be  beneficial  to  focal  plane  array  performance  particularly  in  the 
very  long  wave  infrared  range. 


Oxford,  1979. 

8.  Sher,  A.,  M.A.  Berding,  and  A.-B.  Chen,  J.  Cryst.  Growth,  1989,  98,  27-36. 

9.  Nemirovsky,  A.  and  A.  Ron,  Solid  State  Electronics,  1997,  41,  1881. 

10.  Williams,  G.M.,  R.E.  DeWames,  J.  Bajaj,  and  E.R.  Blazejewski,  J.  Elect.  Mater.,  1993,  22,  931. 

11.  Krishnamurthy,  S.,  T.N.  Casselman,  and  A.  Sher  (unpublished). 

12.  Sher,  A.,  A.  Tsigelman,  E.  Weiss,  and  N.  Mainzer,  J.  Vac.  Sci.  TechnoL,  1990,  A8(2),  1093. 

13.  Eanrie,  J.P.,  J.  Reno,  S.  Sivananthan,  I.K.  Sou,  X.  Chu,  M.  Boukerche,  and  P.S.  Wijewarnasuriya, 
J.  Vac.  Sci.  TechnoL,  1986,  A4(4),  2067. 


